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(N : ABSTRACT 

<' 

^ I We present a medium-resolution spectroscopic survey of late-type giant stars 

^ ■ at mid-Galactic latitudes of (30° < |6| <60°), designed to probe the properties of 

Ci-f this population to distances of ~9 kpc. There are two main goals of this survey. 

O ■ The first goal is to explore the kinematical properties of the thick disk. Since 

■ M giants are generally metal-rich and we have limited contamination from thin 

, disk stars by the latitude selection, most of the stars in the survey are expected 



to be members of the thick disk (<[Fe/H]>~-0.6) with little confusion from the 
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metal-poor (<[Fe/H]>~-1.6) halo. The second goal of the survey is to study the 
properties of the metal-rich component of the nearby halo, a Galactic population 
that is also represented in our sample. 

We report first results for 1799 stars. The distribution of radial velocity (RV) 
as a function of / for these stars shows (1) the expected thick disk population and 
(2) local metal-rich halo stars moving at high speeds relative to the disk, that in 
some cases form distinct sequences in RV-/. High-resolution echelle spectra taken 
for 34 of these "RV outliers" reveal the following patterns across the [Ti/Fe] and 
[Fe/H] plane: sixteen of the stars have abundances reminiscent of the populations 
present in dwarf satellites of the Milky Way; eight have abundances coincident 
with those of the Galactic disk and metal-rich halo; and ten of the stars fall on 
the locus defined by the majority of stars in the halo. The chemical abundance 
trends of the RV outliers suggest that our sample consists predominantly of stars 
accreted from infalling dwarf galaxies, with a more moderate fraction (~20%) 
of stars likely formed in the Galactic disk and subsequently kicked to higher 
eccentricity orbits. These results support scenarios where the stellar halo arises 
from multiple formation mechanisms. We conclude that M giants with large RVs 
can provide particularly fruitful samples to mine for accreted structures and that 
some of the putative velocity sequences may indeed correspond to real physical 
associations resulting from recent accretion events. 

Subject headings: Galaxy: formation - Galaxy: evolution - Galaxy: disk - 
Galaxy: halo - Galaxy: kinematics and dynamics - stars: abundances - solar 
neighborhood 



1. Introduction 

The formation history of the Milky Way is recorded in the present motions and chemical 
abundances of its stars. Ideally, to unravel the Milky Way's history, we would like a catalog 
containing spatial, kinematical, and chemical data for every star in the Galaxy. All-sky 
photometric surveys, such as 2MASS and SDSS, are bringing us closer to this goal: both 



allow for detailed studies of the structure of each Galactic component fe. 




Skrutskie et al. 


2001: Maiewski et al. 2003: Bell et al. 2008: Juric et al. 


2008; 


Ivezic et al. 


2008) and are rich 


sources for follow-up spectroscopic studies (e.g.. 


Majewski et al. 


20o4 


Yannv et al. 2009). 



The present study looks at the spectroscopic properties of a sample of nearby (rf < 9 
kpc) M giant stars at mid-Galactic latitudes of 30° < |6| < 60° selected from the 2MASS 
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catalogue. M giants: are bright stars and, hence, can be easily observed to large distances 
using even small telescopes; can be readily identified on the basis of near-infrared, JHK 
photometry and complete samples can be culled from full-sky catalogs (e.g., 2MASS); are 
limited to more metal-rich populations. These unique properties of M giants, when combined 
with our latitude and magnitude cuts, make our survey particularly useful for exploring the 
structure of the thick disk. Other bright stellar types are not as well-suited for cleanly 
isolating the thick disk, because they span a wider range in metallicity and are distributed 



throughout all parts of the Galaxy. Indeed, iMajewski et al.l ( l2003l ) mapped the distribution 



of M giants in the Galaxy using 2MASS and showed that the majority are found in the bulge 
or disk components of the Galaxy. Moreover, in our survey, contamination from thin disk 
and bulge stars is minimized by the latitude and magnitude limits, and the contribution of 
the stellar halo is expected to be small because these stars have lower metallicities than M 
giants. Nevertheless, we also anticipate some examples of stellar halo stars in our sample: 
Majewski et al. (2003) found a small number of relatively rare, metal-rich M giants in this 
Galactic component. 

Given the properties of our M giant sample, our study has the twin goals of exploring 
the bulk structure of the thick disk and also the nature of the metal-rich nearby stellar halo. 
In this first paper describing our survey, we present the photometric and spectroscopic data 
for our current entire sample of 1799 stars, but focus on the detection and interpretation 
of those stars that do not kinematically conform to the typical behavior of the Milky Way 
thick disk population and are likely to be members of the stellar halo. Our stellar halo 
sample is unique compared to other halo surveys of stars in that it covers an intermediate 
distance range (rf < 9 kpc) and concentrates on a relatively rare halo tracer. For example, 
our survey volume is wider than that of the Hipparcos survey of the Sun's closest neighbors 
(c^iim.Hip ~100 pc - Ferryman et al. 2002) but is a much more local view of the stellar 
halo than the SDSS F-G type turnoff stars (visible to 20-40 kpc - Juric et al. 2008; Bell 
et al. 2008) or the entire 2MASS M giant catalog (visible out to 100 kpc - Majewski et al. 
2003). Our catalog is comparable in scope to the spectroscopic studies of the Century Survey 
Galactic Halo Froject (Brown et al. 2008), RAVE (Steinmetz et al. 2006), and Carollo et al. 
(2007, 2010), although these studies employed different, more commonly used halo tracers 
(BHB stars, stars - both dwarfs and giants - with 9 < / < 12, and MSTO stars and dwarfs 
from the SDSS DR7 calibration sample, respectively). 

In this first analysis of our survey we explore the origin of the M giant population in the 
stellar halo. There are three broad categories of formation scenarios typically postulated for 
these stars: 



1. In-situ-halo stars form at comparable radii to their current locations within the dom- 
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inant dark matter halo progenitor of the Galaxy. For example lEggen et al.l (Il962[ ) 
envisioned the stellar halo could be formed from infalling gas, prior to the forma- 
tion of disk, during an early r nonolithic collapse phase fo r our Galaxy, as seen in the 
hydrodynamical simulation of ISamland &: Gerhard! ( 120031 ) . 



2. Kicked- out- disk stars are formed initially more concentrated towards the center of 
the dominant dark matter halo and are subsequently ejected to ni o re ec centric or- 



bits through minor or major merg ers, as propos e d by jPurce . 



in hydrodynamical simulations by IZolotov et al.l (120091 ) and 



1 et al. (2010) and seen 



McCarthy et al.l tom . 



(Note that there were many earlier studies of this same process that focussed on the 



Quinn et al. 


1993; 


Walker et al. 


1996) 



3. Accreted stars form in separate dark matter halos that later merge with the Galaxy 
(as proposed by Searle & Zinn 1978). 



Note that in prior theoretical work, the first two categories were both simply termed in situ 
to indicate stars formed in the dominant galaxy dark matter halo progenitor (as in, e.g., 
Abadi et al.ll2006l ). Separate terms are introduced here for clarity. 



Many prior studies haye probed these formation mechanisms by looking at the distri- 
bution of stars in different dimensions of phase-space. For example, the importance of an 
accreted population has traditionally been assessed by looking for residual groupings that 
are signatures of stars' original associations: in-situ-halo or kicked- out- disk populations are 
expected to be more smoothly distributed. All-sky photometric surveys have revealed rich 
substructure in the outer halo (e.g., Sgr tidal tails - Majewski et al. 2003; the Orphan stream 
Belqkurov et al. 2007b) that can be best explained by accretion events (e.g.. 



200ll ; iBuUock and Johnstonll2005l ; iBell et al.ll2008l ; Ijohnston et al.ll2008l ; ICooper et al 



Bullock et al. 



20101). 



In contrast, debris in the nearby halo fully phase-mixes on a much shorter timescale, lead- 
ing to the expectation of negligible evidence for accretion identifiable as coherent spatial 
structures (e.g., Johnston et al. 2008; Sharma et al. 2010) and requiring additional di- 
mensions of phase space data to distinguish formation scenarios in this region. Adding the 
dimension of line-of-sight velocities helps: conservation of phase-space density during phase- 
mixing requires that as debris from accretion events becor nes less dense with tim e, it should 
become colder in velocity (Liouville's Theorem; see e.g., iHelmi fc Whitelll999l ) and coher- 
ent structures may be apparent even once stars are smoothly distributed in space. Indeed, 
substructure in velocities is apparent statistically over a large volume in BHB stars in the 
SDSS survey (and shown to be broadly consistent with model stellar halos built within a 
hierarchical cosmology — see Xue et al. 2010, Cooper et al. 2010a) and individual clumps in 
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velocity have also bee n detected more nearby using metal-poor main-sequence turnoff stars 
i n the SEGUE survey ( ISchlaufman et al.ll201ll ). as well as in the mixed populations in RAVE 
( iWilliams et al.ll201l[ ). Overall, group-finding is more effective if even more dimensions of 
phase-space can be measured. For example, Majewski (1992) analyzed proper motions for 
a sample of 250 F-K dwarfs in the direction of the North Galactic Pole that probes out to 
roughly 8 kpc; he measured a mean retrograde rotational velocity for the halo sample and 
detected a more coherent retrograde group of stars at a mean height above the Galactic 
plane of Z ~ 4.6 kpc - findings suggestive of an accreted halo population. Subsequently, 
Majewski, Munn, & Hawley (1994; 1996) obtained spectra for a subsample of stars from the 
Majewski (1992) survey and found a significant amount of dumpiness in their halo sample. 
Once all six phase-space dimensions are known, conserved (or nearly conserved) quantities 
(such as energy, angular momentum or orbital frequencies) can be calculated that can link 
stars from a common progenitor in the volume even if they are no t at the same orbital 



phase, as has beeri done successfully for seve r al nearby surveys (see iHelmi fc White! Il999 
Helmi fc de Zeeuwiboool : iMorrison et"aDl2009l : lOomez et allboioh . 



While the findings in the previous paragraph point to a significant fraction, and possibly 
the majority of the stellar halo being accreted, some contributions from in-situ-halo and 
kicked- out- disk stars are still possible. Structural, orbital and/or metallicity trends in the 
stellar halo with radius — as well as transitions in those trends — have historically been taken 



and b lue-horizontal -branch stars, se e 



Hartwick 


1987; 


Zinn 


1993, 


1996: 



19971 . respectively). IChiba fc BeersI (|2000[ ) analyzed a kinematically unbiased sample of 1203 
stars with [Fe/H]<-0.6 in the inner halo (within 4 kpc of the Sun) and detected a negative 
rotational gradient as a function of height above the Galactic plane for the more metal- 
rich stars in their sa r nple — as seen for the in-situ-halo stars formed in the simulations of 
Samland fc GerhardI ( l2003l ). They also confirm the existence of the streams detected by 
Helmi et al.l (119991 ) — indicative of some presence of an accreted population in this region. 
Carollo et al. (2007, 2010) studied the kinematics and metallicities of SDSS calibration stars 
for a larger volume (20 kpc) and found evidence for two populations: one of metal-rich stars 
on only mildly eccentric orbits (which they dubbed " inner halo" ) and a second of metal-poor 
stars on more eccentric orbits (which they dubbed "outer halo"). Similarly, iDeason et al. 



( I2OIII ) analyzed BHB stars in SDSS (a sample that probes the outer halo to 40 kpc) and 
find a net retrograde rotation for metal-poor stars and a net prograde rotation for metal-rich 
stars. Such mult i- component halos, with distinct formation mechanisms for each component 
emerge naturally in the hydrodynamical simulations, with the inner halo (within 10-15 kpc) 
coming predominantly from stars (either in-situ-halo or kicked- out- disk) formed within the 
main Galaxy dark matter halo progenitor and the outer halo (dominant beyond 15-20 kpc) 
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from mergers and accretion of stars by the Galactic dark matter halo (lAbadi et al.l 12006 
Zolotov et ahlboogi IPont et al.lboill : iMcCarthv et al.ll2012[ ). However, whether the distinct 
components observed in the stellar halo are indeed due to distinct formation mechanisms, 
and not merely different accretion events, has yet to be proved. 

In our own survey, the sample of stars is sufficiently far that distance and proper motion 
errors are too large to estimate their orbital properties accurately. However, chemical abun- 
dances can be derived for the high-resolution spectra, thus providing an alternate avenue 
for exploring the origins of the stars in our survey. A star is branded at its birth by its 
chemical abundance — a signature that is generally conserved throughout its lifetime (like 
orbital properties) and cannot be diluted by orbital phase-mixing. Moreover, stars deriving 
from a common origin should have similarities in their chemical abundance patterns, trends 
that are directly correlated to the details of their enrichment history. Hence we can hope to 
"chemically tag" stars as members of the different populations via their abundance patterns 
(similar in spirit to Freeman & Bland-Hawthorn's 2002 proposal for the reconstructing of 
ancient star-clusters in the stellar disk). The potential power of chemical tagging has already 
been demonstrated empirically in observations: in dwarf galaxies, stars tend to have lower 
[a/Fe] at a given [Fe/H] than stars in the bulk of the Milky Way's stellar halo (Smecker-Hane 
& McWilliam 2002; Shetrone et al. 2003; Tolstoy et al. 2004; Geisler et al. 2005; Monaco et 
al. 2007; Chou et al. 2010a) and similar patterns have been seen in stars in str uctures such 



as th e Monoceros Ring (Chou et al. 2010b) and Triangulum-Andromeda Cloud (IChou et al. 



20111 ). which lends support to the interp retation of such features as disrupted dwarf galaxies . 
In a similar manner, a series of papers ( iNissen &: Schusterll2010l . 1201 ll : ISchuster et al.ll2012[ ) 
have separated a sample 94 nearby (within ~335 pc of the Sun), metal rich stars (-1.6 < 
[Fe/H] < -0.4) into "a-rich" and "a-poor" groups and shown systematic differences in the 
abundances, ages and orbital properties of stars in these two groups that are suggestive of 
kicked- out- disk and accreted origins respectively. 

Figure [1] illustrates how this approach can be applied to our own sample by showing 
where populations of in-situ-halo, kicked- out- disk and accreted stais might lie in the [a/Fe]- 
[Fe/H] plane given our expectations of systematic differences in each population's history. 
For example, in the hydrodynamical simulations, the kicked- out- disk population forms in 
the deep potential of the primary Galactic progenitor prior to the last major merger event 
( IZolotov et al.ll2009l : [McCarthy et al.ll2012l ). and hence their abunda nces may be expecte d to 
overlap with the older, thick disk populations (as first pointed out by lZolotov et al.ll2010l ). In 
contrast, stars in accreted populations (shown in Figure [T]as the green ellipse) form in smaller 
potential wells (by definition) over the entire history of the Milky Way (e.g.. Bullock & 
Johnston 2005). Star formation and metal enrichment are less efficient in these smaller dark 
matter halos, allowing enrichment by Type la SNe — and evolution of [a /Fe] to low values — 
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befor e the mean metallicity has risen as high as the disk ( iRobertson et al.ll2005l : iFont et aL 



20061 ). Note that stars formed via all three mechanisms {in-situ-halo, kicked- out- disk or 
accreted) could have low metallicites and high [a/Fe], as seen in the bulk of the stellar halo 
(shown as the blue ellipse) and thus we do not anticipate being able to definitively separate 
out all three populations in this region of abundance space. Nevertheless, these expectations 
lead us to conclude that the fraction of stars in a metal-poor and a-poor population is 
an indicator of the importance of late accretion, while the fraction of stars with disk-like 
abundances but moving at large speed relative to the Sun is indicative of the contribution 
of a kicked- out- disk population. 

An additional advantage of focusing on abundance space as a probe of the origin of the M 
giants in the stellar halo is that we only need to look for the expected systematic differences 
in the chemical composition of these populations [in-situ-halo, kicked-out disk, and accreted) 
overall rather than (for example) search for kinematical groupings from individual accretion 
events. Hence, we can explore formation scenarios with much smaller samples of stars than 
by looking at dynamics alone. Motivated by the promise of chemical abundances, in addition 
to radial velocities for all program stars, we obtained high-resolution echelle spectra for 34 
stars that were selected for follow-up based upon their high, halo-like speeds relative to the 
disk, as well as five random thick disk red giants (cho sen based on RVs simila r to the bulk 



thick disk trend) and four M giant calibration stars (jSmith &: LambertI Il985l : ISmith et al. 



2000). 



This paper is organized in the following way. In §21 we describe target selection, ob- 
servations and data reduction for the medium-resolution spectroscopy program of bright M 
giants. The RVs and the identification of several distinct groups of RV outliers - stars that 
have high speeds relative to the bulk thick disk trend - are presented in ^|3l Details of the 
high-resolution follow-up spectroscopy program are given in §H Lastly, we give a summary 
of the results and disc uss them in the context of prior studies in §3 In a companion paper 



(I Johnston et al.ll2012l ). we go on to develop a more complete understanding of the nature 
and implications of the chemo-dynamical groups in our survey by generating and analyzing 
synthetic observations of simulated stellar halos. 



2. Program Stars 
2.1. Defining the Sample 



M giants can be distinguished fr om M dwarfs by their J — H and H — K colors 
(iBessell &: BrettI Il988l : ICarpenterl 1200 ll ). M giants also dominate M dwarfs in catalogs of 
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late-type stars to roughly K ~14. These facts were used by iMajewski et al.l (120031 ) to se- 
lect M giants from 2MASS and map t he streams of M gi a nts from the Sagi ttarius (Sgr) 
dwarf spheroidal (dSph) galaxy and by iRocha-Pinto et al.l (l2003l . |2004| . |2006| ) to map the 
Triangulum- Andromeda star cloud and to track the Monoceros/GASS/Argo feature. Our 
survey sam ple spans {J—Ks)n c olors from 0.75 < {J—Ks)o < 1.24; this is similar to the range 
studied by iGirard et al.l (120061 ) in their study of the thick disk using red giants, although 
most (96%) of our red giants have (J — -^'5)0 > 0.85 (as in Majewski et al. 2003) to ensure 
a clean sample of M giants. The magnitude range of our sample is 4.3 < (-^^5)0 < 12.0, with 
a mean and dispersion of 7.5±1.3 (the majority, 1631 of 1799 have 5.0 < (-^'5)0 < 9.0). The 
magnitudes were dereddened using the maps from Schlegel et al. (1998), and we removed 
from our sample any stars with E{B ~V) > 0.555. 

Because we intended to sample primarily the thick disk, a constraint in Galactic lati- 
tude of 30° < |6| < 60° was applied to our M giant catalog: the lower limit in h is meant 
to avoid excess contamination from the thin disk and the upper limit is intended to further 
reduce inclusion of halo field stars. Our nominal photometric sample contained approxi- 
mately 12000 stars, and to date we have spectroscopically observed 1799, or roughly 15% 
of these. The stars observed were selected randomly from the nominal sample and cover 
nearly all Galactic longitudes. The bulk of the observing was done at the Fan Mountain 
Observatory (FMO), located in Virginia, so there are gaps in coverage corresponding to the 
Southern Hemisphere. Of the 1799 stars, 149 were observed at Cerro-Tololo Inter-American 
Observatory (CTIO) as part of a related program. Figure [2] shows the spatial distribution 
in an Aitoff projection of the 1799 program stars in Galactic coordinat es. The nominal M 
giant catalog was matched to the UCAC2 catalog (IZacharias et al.ll2004l ) and the gap in the 
northern Galactic hemisphere is due to the upper limit oi 5 = +52° in the UCAC2 catalog at 
the time of the survey inception. However, the very small amplitudes of the UCAC2 proper 
motions for the majority of the program stars typically result in unreasonably large relative 
errors (larger than the space motions) in any derived kinematical parameters, so the proper 
motions are not in cluded in this ^york. L imitations using the UCAC2 proper motions are 
further explored in IMajewski et al.l (120061 ). 



2.2. FMO and CTIO Spectroscopic Observations and Reductions 



Spectra were collected at the University of Virginia's FMO using the Fan Observatory 
Bench Optical Spectrograph (FOBOS; see Crane et al. 2005) on the 1-m astrometric reflec- 
tor. FOBOS is a fiber-fed opt ical spectrograph that was designed for moderate resolution 
{R ~ 1500-3000) spectroscopy (ICrane et al.ll2005l ). FOBOS uses a grating with 1200 grooves 
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mm~^. The estimated resolution of the spectra is AA ~ 4 A. Our observing program began 
on UT 2005 February 25. The spectrograph is optimized for use over the region 4000-6700 
A but is hmited by the camera, which has a SITe 2048x2048 CCD detector, which at a 
hnear dispersion of 1.0 A pixel"^ samples ~2000 A of that range (4000-6000 A). On all but 
one night, the detector setting used had a read noise level of 4.5 e- and a gain of 6.1 e- 
per ADU. For most stars, three spectra are taken and summed in 2-D after the CCD frame 
preprocessing is completed. This combination of three images facilitates the elimination of 
cosmic rays. The minimum S/N to achieve the best possible RV precision (a few km s^^) 
was found to be ~20; this typically translated to total exposure times of 450-900 seconds for 
the M giants observed with FOBOS, which have magnitudes in the range 4.3 < (-^'5)0 < 9.7, 
with a mean value and dispersion of 7.3±1.0. Several radial velocity standard stars from 
the HD catalog were also observed each night; these are used for cross-correlation templates 
when determining the radial velocities of the target stars. Standard stars were chosen to be 
of a similar spectral type as the program stars to minimize systematic offsets in the derived 
RVs. Several sets of bias frames were taken throughout each night. To remove pixel-to-pixel 
variations in the frames, a "milky flat" is created by illuminating an opal diffusing glass with 
a quartz-tungsten-halogen lamp. The object and comparison frames are fiat-fielded (using 
the milky fiat), trimmed, and bias subtracted using the task ccdproc in Comparison 
spectra were taken using neon, argon and xenon lamps for calibrating wavelengths against 
the laboratory values for lines from these elements. The spectra are extracted from the 2-D 
images and converted to 1-D spectra and wavelength calibrated using the IRAF tasks apall 
and identify. 

To obtain radial velocities, the extracted, fiat-fielded, wavelength-calibrated spectra are 
cross-correlated with the standard star spectra using code developed by W. Kunkel (described 
in detail in Majewski et al. 2004). The cross-correlation code is run for each night of data 
and each program star is cross-correlated to all of the standard stars observed on that night. 
The radial velocity reported for a program star is the average of the radial velocities from 
cross-correlation with multiple standards taken that night (but standard star spectra that 
produce poor cross-correlations against the others are removed from the average, iteratively). 

The spectrum for a typical program star observed at Fan Mountain Observatory is 
shown in Figure El The dominant features in M giants in the spectral band we are studying 
are the Mg b triplet, at 5167, 5173, and 5184 A, and the Na D doublet at 5889 and 5896 
A. A number of strong Fe, Cr, and Ti lines/blends are also present. The redder stars in 
our sample (J — Ks > 1.1) show very strong titanium oxide (TiO) bands in their spectra; 



^IRAF (Image Reduction and Analysis Facility) is distributed by NO AO, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under contract with the NSF. 
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these M giants are cool enough (Teg < 3560 K) that the TiO bonds in the star are not 
dissociated. Figure H] shows the spectra for several stars observed with FOBOS covering a 
range of (J — Ks) colors; the appearance of the TiO bands is apparent as the giants become 
redder and cooler. 

An additional 149 M giants that fit our survey criteria were observed over UT 2004 
October 8 - 11 at CTIO using the Cassegrain spectrograph on the 1.5-m telescope. The 
detector was a Loral IK (1200x800 pixels) CCD with a read noise of 6.5 e-; the gain was 
set to 1.42 e- per ADU. A grating with 831 grooves mm~^ was used, with a resolution of 
A A ~ 3.1 A in the first order. Helium and argon lamps were used to take comparison 
frames for each target at the same telescope position to account for flexure variations. Ten 
quartz frames were taken each night and combined and normalized. The spectral range is 
7650-8900 A. The dominant feature in this region is the Ca II triplet at 8498, 8542, and 8662 
A. The CTIO M giants have magnitudes in the range 6.3 < (-^'5)0 < 12.0, with a mean and 
dispersion of 9.8±1.2. The reduction procedures for the CTIO program giants are similar to 
those used for the FMO reductions (the same cross-correlation code was used to determine 
the radial velocities). 

For the FMO sample, the difference in the derived (from cross- correlation against each 
other) and published RVs for the standard stars is typically ±1-5 km s~^, with no systematic 
offset in either direction. A total of 102 program stars were observed multiple times at FMO, 
to test the stability of FOBOS and to gauge the S/N threshold for obtaining reliable RVs. 
The mean of the absolute value of the differences in the FMO RVs is 8.5 km s~^. The RVs 
are fairly stable even at low S/N; stars with S/N below 20 have a mean absolute value in the 
difference of their RVs of 10.5 km s~^. For the CTIO sample, the difference in the derived 
and published RVs for the standard stars is slightly higher, with variations ranging from 
±1-10 km s~^; as with the FMO standards, no systematic offset in either direction is seen. 
Repeat observations were also taken for 16 stars at CTIO. The mean of the absolute value 
of the differences in RVs for the CTIO repeat observations is 13.9 km s~^. 



3. Radial Velocity Distribution 



Panel (a) of Figure [5] shows the heliocentric radial velocity, function of Galac- 

tic longitude, /. In panel (b), these velocities have been translated to the GSR (Galac- 
tic Standard of Rest frame — i.e., centered on the Sun but at rest with respect to the 
Galactic Center), where we adopt the values Oo=236 km s^^ for the speed of a c losed 
orbit at the position of the Sun relati ve to the Galactic ce nter (IBovy et al.l l2009l ) and 
([/0,1/0,W^0) = (11-1, 12.24,7.25) km s'^ JSchonrich et al.l boioh for the motion of the Sun 
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with respect to this orbit. General trends in these panels can be understood by assuming 
that stars in the Galactic disk move on nearly circular orbits around the Galactic center. 
For a flat Milky Way rotation curve near the Sun with circular speed Gq, the predicted 
line-of-sight velocity with respect to the GSR for a star on a circular orbit is given by: 



Here, R is the Galactocentric radius of the star and Rq is the solar Galactocentric radius. 
The expected sinusoidal trend for stars moving on disk-like orbits around the Galactic center 
is seen for the bulk of the stars. (Note that the heliocentric velocities in panel (a) still show 
some sinusoidal trend — a reflection of the mean motion of the Sun with respect to the M 
giant population. )@ 

At fixed Z, equation ([T]) shows that fGSR,circ is lower for stars on circular orbits observed 
at high h than for stars observed at low h due to the cos 6 projection of their motions. In 
principle, therefore, one may have a good approximation to the rotational component of 
a disk star's velocity at a latitude h by deprojecting the observed velocity with a division 
by cos 6. This scaling tends to accentuate differences between stars having "disk-like" (i.e., 
circular) motion from stars having non-disk-like motions because the deprojection will tend 
to tighten the coherence of disk stars but separate outliers more — as shown in panel (c) of 
Figure |5l 

While the majority of the stars in Figure [5] appear to be members of the thick disk, 
there are a significant number of stars with high velocities relative to the main trend (we 
refer to these stars as "RV outliers"). The origin of these RV outliers — in-situ-halo, kicked- 
out-disk or accreted — is unclear from this plane alone. If the RV outliers are in-situ-halo 
or kicked- out- disk stars, we would expect to see random RVs as a function of I. However, 
stars in several longitude ranges in our sample show suggestive coherence in their RVs. Such 
coherence is expected for a stream of stars passing through the Solar neighborhood (see 
Johnston et al. 2012 for further details), but whether these structures are real is hard to 
assess given the small number of stars. 

In the next section we present high-resolution follow-up spectroscopy of RV outliers to 
further explore what fraction of outliers can be attributed to each of the halo formation 
mechanisms. 




(1) 



^The full medium-resolution radial velocity catalog is available upon request. 
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4. Chemical Abundances 

4.1. Data Collection and Reductions 

We have assigned the RV outliers to groups by eye, based upon apparent correlated 
sequences in RV-Z space as indicated in Panel (d) of Figure [51 To test the chemical properties 
of the stars in some of the RV outlier groupings, high-resolution echelle spectra were collected 
using the Astrophysical Research Consortium Echelle Spectrograph (ARCES) on the 3.5-m 
telescope at Apache Point Observatory on UT 2009 March 30 and UT 2009 April 2 and 
using the CCD echelle spectrograph (ECHLR) on the 4-m Mayall telescope at the Kitt Peak 
National Observatory over UT 2010 February 26 - March 2. 

The ARCES uses a 2048x2048 pixel SITe CCD and has a resolution of i?=31,500; the 
CCD has a gain of 3.8 e- per ADU and a readout noise level of about 7 e-. Two sets of fiat 
fields were taken to account for the strong gradient in response across the ARCES orders: 
one long set with a blue filter inserted and another short set with no filter (these are the red 
calibration frames). The blue and red frames were combined to create a master quartz flat. 
For wavelength calibration, thorium-argon (ThAr) lamp frames were taken. 

Reduction of the ARCES data was carried out using various IRAF tasks in the echelle 
package. All images were overscan-corrected and trimmed using ccdproc. The echelle orders 
were located and the trace defined for the spectra with apall. The task ecidentify was 
used to identify lines in a ThAr lamp spectrum. To minimize the effects of aliasing, the 
spectra were resampled along the dispersion axis using the magnify task. Scattered light 
was removed from the program star frames using the apscatter task, and the relevant echelle 
orders were then extracted. The program spectra were divided by the extracted master 
quartz fiat, and the dispersion correction defined using the ThAr lamps was then applied to 
convert the pixel scale to a wavelength scale. As a final step in the reduction process, the 
spectra were continuum normalized. 

ECHLR spectra were collected at KPNO using the 2048x2048 pixel T2KB CCD on the 
4-m Mayall. The gain setting for T2KB was 1.9 e- per ADU with read noise approximately 4 
e-. The ECHLR has a resolution of i?=35,000. ThAr lamps were used to collect wavelength 
calibration frames, and the short-exposure quartz lamp was used for obtaining the fiat-field 
images. The spectral reduction procedures for the ECHLR data is similar to that carried 
out for the ARCES data. 

The photometric properties and the observational details of the stars observed at APO 
and KPNO are listed in Table [H The IDs and photometry of the stars come from the 2MASS 
point source catalog, where the IDs are the 2MASS RA/Dec (J2000.0) coordinates. Radial 
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velocity standards, taken from the Astronomical Almanac (2008, 2009), were observed all 
nights at APO and KPNO. Radial velocities from the echelle data were found using the fxcor 
task and are reported in Table [Hand are compared with the medium- resolution values (high- 
resolution/low-resolution, denoted as Whei,h/'i'hei,m5 respectively). Cross-correlation between 
RV standards gives errors on the order of 0.5 km s~^ for the high-resolution RVs. The 
mean of the absolute value in the differences between the medium-resolution RVs with the 
high-resolution RVs is 6.2 km s~^. Based on the RV data in Tabled! there is apparently a 
systematic bias towards higher measured RVs for the medium-resolution spectra, such that 
^^hei,h — ^^hei,m=-5.0 km s~^. Overall, considering the random errors in the medium-resolution 
RVs and the comparison of the medium/high resolution values, we place the uncertainty 
level for the medium-resolution RVs at 5-10 km s~^. 



4.2. Chemical Abundances 



4-2.1. Derivations 



The APO and KPNO instrument set-ups give the best S/N per pixel for the region 
around 7400 A. This particular region of the spectrum was selected due to its relative absenc e 
of molecular bands, offering a good window for spectral analysis ( ISmith fc Lambertlll990l ). 
The procedures used to derive the metallicities are similar to those used by Chou et al. 
(2007) in their study of M giants in the Sgr tidal tails. Equivalent widths for Fe I lines were 
used to determine the atmospheric parameters Tgfj, log g, and [Fe/H]. The equivalent widths 
were measured manually using the IRAF splot task for 11 Fe I lines in the range 7443 to 
7583 A. 

The measured equivalent widths (EW) comprise a line list used as input to the MOOG 
local thermodynamic equilibrium code (Sneden 1973). The wavelengths of the 11 Fe I lines 
and the corresponding EW measured for each star are listed in Table |21 In some cases, 
a line could not be measured due to a cosmic ray falling on the same pixel as an Fe line; 
these lines were removed from the list for that spectrum. We adopt the same values for 
the excitation potentials (x) and oscillator strengths [gf) for these lines as those used by 
in Chou et al. (2007). Along with the line list, stellar atmosphere models from the Kurucz 
ATLAS 9 grids (1993) are used. Each stellar atmosphere model corresponds to a particular 
(Teff, log g, [Fe/H]). We start with an initial guess for (Tpff, log g, [Fe/H])o. The starting 
value for Tes is informed by the iHoudashelt et al.l ( l2000l ) color-temperature relations for M 
giants, usin g the dereddened 2MASS J — Ks colors converted to the CTIO system and the 
relations of ICarpenterl (120011 ). An initial guess for log g is taken from the 5 Gyr Zq isochrone 
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from the Padova evolutionary track database (Marigo et al. 200qj). The surface gravity g 
of a star is related to its mass and size (luminosity). Once a star exhausts the hydrogen 
supply in its core, it will first brighten and move up the red giant branch (RGB); later, after 
core He exhaustion, the star eventually ascends the asymptotic giant branch (AGB). The 
difference in log g between a star on the RGB versus the AGB for cool giants, however, is 
quite small over a wide range of ages. Figure |6l which shows the Teff-log g plane for stars of 
age 3 Gyr, 5 Gyr, and 10 Gyr for a star of solar metallicity (^©=0.019), demonstrates that 
the variation in log g between the RGB and AGB is on the order of 0.15 dex for stars in the 
range of effective temperatures for our sample, which is roughly 3600-4000 K. The iterative 
scheme — run MOOG, update (Tcs, log g, [Fe/H]) based on MOOG results, select updated 
stellar atmosphere model from the ATLAS 9 grids — is continued until the slope of the Fe 
abundances as a function of log(x/A) is less than ±0.005 dex. The derived parameters for 
all program stars are listed in Table [31 

In addition to [Fe/H], the ratio [Ti/Fe] was also derived. Once the atmospheric param- 
eters for a star were determined using MOOG, these were fixed and used to measure A(Ti). 
Three Ti I lines were used to derive A{Ti): 7474.940, 7489.572, and 7496.120 A. For the first 
line, the oscillator strength was derived by fitting the EW of the line to the solar spectrum. 
The oscillator strengths for the latter two Ti I lines are taken from Chou et al. 2007 and 
were derived using the same technique of fitting the EW of the line to the solar spectrum. 
A s a check of our methodo logies , four calibra t ion s tars from the red giant spectral studies 



of ISmith &: LambertI (Il985[ ) and ISmith et al.l (120001 ) (collectively referred to as S&L) were 



observed with both the ARCES and ECHLR; a comparison between our [Fe/H] and [Ti/Fe] 
values for the calibration stars is given in Table |H There is a notable discrepancy in the 
derived [Fe/H] values for u Vir. Typical uncertainties in the values derived for [Fe/H] in both 
studies are ~ ±0.15 dex using this sample of Fel lines. The mean difference in [Fe/H], and its 
associated standard deviation, between these two studies, in the sense of (this study - S&L) 
= -0.04 ± 0.16. This comparison indicates that the two studies are on the same abundance 
scale, with only a small mean offset, and a scatter typical of the derived uncertainties in 
[Fe/H]. The values are in good agreement and the discrepancies are due to differences in the 
stellar atmosphere models used. 



see http://pleiadi.pd.astro.it/ 
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4-2.2. Distances 

High-resolution spectra allow us to compute distances. These distances are helpful in 
this paper to obtain some understanding of the size of the volume our M giant sample 
probes around the Sun as well as their height above the Galactic plane. The distances to 



(Marieo et al. 


2008; 


Girardi et al. 


2010) 



derived [Fe/H] and assuming stars of age 5 Gyr; the appropriate Mk is then selected based 
on the derived Tes- Using an age of 5 Gyr means that a star evolving along the RGB or AGB 
would have had an initial mass of ~ 1.1-1.3 Mq; such a typical mass would not be largely 
different from a much younger population, such as 2 Gyr where Mrgb/agb ~ I.TMq, or a 
much older population with an age of 10 Gyr and Mrgb/agb ~ I.OM0. As discussed above 
(see Figure [6]), the stellar atmospheric parameters for a red giant on the RGB do not vary 
significantly from those for a red giant on the AGB. 

The distances to all program stars are listed in Table O the derived distances of the 
standards agree well with their Hipparcos distances (10% to 14% accuracy) when using 
isochrones in the same way as used for the program stars. Possible sources of systematic 
error in the derived distances include: (1) an incorrect assumption of the star's age and (2) 
scatter in the isochrone for later evolutionary stages. As a check on the error in distance from 
the isochrone method, was computed for ages of 2, 5, and 8 Gyr for stars having different 
values of Tgfj and [Fe/H]. Using estimated uncertainties in Tgfj of ~100 K and ~0.10 dex in 
[Fe/H], a conservative distance uncertainty (from the variation in Mk with age, Tgfj, and 
[Fe/H]) is ~25%, with the scatter defining a distribution of distances (such as a) expected 
to be somewhat smaller than this. 

We note that a handful of the more metal-poor stars cannot be fit to an 8 Gyr isochrones, 
as their values of T^s are cooler than the RGB tip. It is possible that this subset of stars are 
cool AGB stars, with the low-metallicity, low-mass stellar models not modeling the end of 
the AGB well. 

Based on the errors in the derived distances for the standards and the sources mentioned 
above, we assign an approximate accuracy of 20% to the distances for the program M giants. 
The distances for the stars with high-resolution spectra confirm that the survey sample 
probes out to roughly 9 kpc from the Sun. 
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4-2.3. Results & Interpretation 

Figure [7] shows [Fe/H] versus [Ti/Fe] for the program M giants having high-resolution 
spectra with points color-coded according to apparent RV outlier groupings (the same color 
scheme as in Figure |5]); two stars with high errors in their derived [Ti/Fe] are not included. 
Galactic disk and halo stars taken from various sources in the literature (Fulbright 2002; 
Reddy et al. 2003; Brewer & Carney 2006) are indicated by open black symbols; the stars 
shown in yellow are presumed thick disk stars with high-resolution spectra from this survey 
classified as such based upon RVs that fall along the bulk disk trend (Figure [5]); as may be 
seen in Figure [7] these stars all fall along the expected main Galactic [Ti/Fe]- [Fe/H] trend. 
The solid line in Figure [7] is the linear [a/Fe] trend for "halo-like" stars while the dotted line 
is the [Fe/H] at which [a/Fe] begins to transition to a decline toward lower values. 

The locations of our RV outliers in Figure [7] can be interpreted in terms of stellar halo 
formation scenarios: 

• Ten of our RV outliers fall along the main halo trend (to the left of the dotted line 
and highlighted by diamonds) but could be in-situ-halo, kicked- out- disk or accreted. 
Because of the expected chemical overlap of stars from these formation scenarios, the 
origin of stars in this region is ambiguous. 

• Eight RV outliers fall along the high-metallicity disk trend (to the right of the dotted 
line) and could be part of the kicked- out- disk population (highlighted by outlining 
squares in Figure [7]), based on the combination of disk-like chemistry but extreme 
kinematics. 

• Sixteen RV outliers fall outside both the disk or halo chemical trends and are likely to 
be part of the accreted population (highlighted by outlining triangles in Figure [7]). 

The population assignments for the program M giants are given in the final column of Table 
1; the column heading of 'o' denotes origin, where 'is' refers to in situ, 'ko' refers to kicked- 
out-disk, and 'a' refers to the accreted population. The five stars with an origin of 'd' are red 
giant thick disk calibration stars, selected based on the fact that their medium-resolution 
RVs fall along the main thick disk RV trend in Figure |5l 

The fraction of M giants in our RV-outlier sample that are potentially members of the 
kicked-out-disk population is small but not negligible (8 out of 34). Figure [H] summarizes 
our measurements of the radial velocity and distance from the Galactic plane for our high- 
resolution sample, emphasizing that these stars have a low probability of being genuine disk 
members both because of their large RVs and because they typically sit at several thick-disk 
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scale-heights (assumed to be 750 pc and shown as a dotted line; ISiegel et al.ll2002l ) away from 
the plane. We can derive a crude upper-limit to the contamination we expect from genuine 
disk stars in our RV-outlier sample by ignoring the decrease in disk density at high-|Z| and 
simply asking what fraction of thick disk stars could be moving at high enough speeds to be 
in our outlier sample. If the 1799 M giants with medium resolution spectra were all drawn 
from a gaussian, with velocity dispersion 50 km s~^ (for a thick disk population), we would 
expect 1.2%, or 22 stars, to have RVs greater than 125 km s~^. In fact, there are 167 stars 
in our medium-resolution sample with such high RVs, suggesting that the majority of these 
are in fact members of another population (i.e., the stellar halo), and that 22/167=13.2% is 
an upper limit to the fraction of stars with RVs greater than 125 km/s that are Galactic disk 
members. Twenty-six out of the 34 high- resolution spectra have RVs (|f/iei|) greater than 125 
km s~^ (dashed lines in Figure [HD, and 4 of these (or 18%) have abundances consistent with 
the Galactic disk, while another 2 with disk-like abundances (stars with IDs 0807206+435418 
and 1313176-1642203) lie within a few km of our velocity cut. 

Overall, we conclude that RV outliers in our M giant sample appear to be dominated 
by accreted stars with a small but significant contribution from potentially kicked- out- disk 
stars. 

It remains to be determined whether any groupings identified by eye in the RV-/ plane 
are actually real. However, some of our outlying stars have RVs and positions that coincide 
with re cent detections of sub structure in the direction of Aquarius at high fhei in the RAVE 
survey (jWilliams et al.ll201ll). In Figurep, we c ompare the location of our stars in RV-/ space 
with the region in which IWilliams et al.l (120111 ) detected the Aquarius stream. In panel (a) 
of Figure [9], f hei as a function of / is shown for only the program stars with negative Galactic 
latitudes (to mimic the RAVE sample selection); panel (b) shows all program stars, with 
the stars having high-resolution spectra plotted as the colored circles. The box extending 
from 30° < / < 75° at negative RVs indicates the r egion in RV-/ space o f the detection of the 
main Aquarius stream (studied most extensively IWilliams et al.ll201l[ ). Note a coincident 
apparent sequence in our southern-hemisphere M giant population in this Aquarius stream 
region of parameter space. Based on a dynamical model of a satellite dissolving in the 
Galaxy's potential, Williams et al. (2011) speculate that the Aquarius stream may be the 
remnants of tidal disruption due to a recent merger. Their interpretation is consistent with 
the notion that some fraction of the giants in our survey with large RVs relative to the disk 
represent recently accreted populations. 
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5. Summary, Discussion and Conclusion 
5.1. Key Results 

In this study, we have analyzed the spectroscopic properties of a sample of M giants 
that are dominated by members of the Milky Way's thick disk and nearby halo components. 
From the RV distribution, we identified stars with RVs that lie outside those expected for 
typical thick disk stars at the same locations. These RV outliers are found to show some 
degree of spatial and kinematical coherence. We suppose that some of this coherence could 
be the signature of substructure (e.g., tidal tails) from accreted satellites. 

To test our interpretation of the RV outliers, we looked at the chemical abundance 
patterns of 34 of these stars. We used the locations of the stars in the [Ti/Fe] versus [Fe/H] 
plane to attempt to assign them to one of three expected populations of halo stars - those 
formed in situ, those that were accreted, and those that were kicked out of the disk. This 
cannot be done unambiguously in all cases because of overlaps in the chemical signatures. 
However, the [Ti/Fe] abundances for sixteen of the RV outliers are systematically below 
the main halo trend, and similar to the abundances seen in stars from Milky Way dSphs. 
They are consistent with stars from MW dSph galaxies that have been accreted by the Milky 
Way (Shetrone et al. 2003; Tolstoy et al. 2004; Monaco et al. 2005) and inconsistent 
with expectations for abundance patterns from pure in-situ-halo or kicked- out- disk growth. 
Moreover, in several cases, stars that were close kinematically also showed some chemical 
coherence, further supporting our accretion hypothesis. Another eight of the stars selected 
for high-resolution spectroscopic follow-up track the abundance trends set by the disk stars, 
even though they have halo-like kinematics. These are indicative of a population formed in 
the disk and subsequently kicked-out, although we estimate some contamination by genuine 
disk stars even at these high RVs and distances from the Galactic disk. The remainder 
of the stars have high [Ti/Fe] at low [Fe/H] and could be part of any of the in-situ-halo, 
kicked- out- disk or the accreted populations. 



5.2. Our results in the context of other studies 



While it has had a long history (jHartwicklll987l : IZinnlll993l . Il996l : ISommer-Larsen et al. 



19971 ). the discussion of possible multiple mechanisms for the formation of the stellar halo 
has been revitalized in recent years both by the ad vent of large photometri c data sets 
with follow-up low-resolution spectroscopic work (e.g.. ICaroUo et al.l 120071 . l2010l ). and more 
modest, but ever-expanding samples of halo stars with high-resolution spectroscopy (e.g., 
Nissen fc Schusteiil2010l . 1201 ll : ISchuster et al.ll2012l ). In several cases, the studies point out 
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classes of halo populations with distinct properties that are argued to match expectations 
for the properties of populations with distinct formation histories seen in hydrodynamical 
simulations of galaxy formation: the models general predict an inner halo (within ~20-30 
kpc) dominated by metal-rich stars formed within the main galactic progenitor, and an outer 



halo dominated by lo wer-metallicity, accreted stars (lAbadi et al.l l2006l : IZolotov et al.l 12009 



McCarthy eraPboid l However, comparisons up to this point are necessarily inconclusive 
since the results of the simulations themselves are dependent on the (hard to model!) details 
of when, how and where stars form in main galactic progenitor. 



For example, ICaroUo et al.l (120071 . |2010| ) found evidence for two populations in their 
sample of 10123 thousands of nearby (within 4 kpc) SDSS calibration stars: one of metal- 
rich stars (peaking around [Fe/H]= -1.6) on only mildly eccentric orbits and a second of 
metal-poor stars on more eccentric or bits These observati ons are very reminiscent of the hy- 
drodynamical simulation results, and ICaroUo et al.l (120101 ) interpret their metal-rich stars as 
an inner, in-situ/ kicked-out-disk population, and their metal-poor stars as an outer, accreted 
population. However, a transition in the orbital and metallicity properties in populations is 
not necessarily inconsistent with a purely accreted stellar halo. Moreover, their interpreta- 
tion is at odds with our finding of a large fraction of clearly acc reted stars in our M giant 
sample of even higher metallicity stars ([Fe/H]> —1.2) than the ICaroUo et al.l (120071 . 12010) 
metal-rich population. In fact, we might expect the M giants to be biased towards finding 
in- situ /kicked- out- disk stars. However, the selection of RV outliers would bias our sample 
towards the high velocity tails of all populations — since accreted stars are expected typically 
to be on higher energy and higher eccentricity orbits than the other stars, this could explain 
why our sample is particularly sensitive to the accreted population. Further work is needed 
to definitively show which of these two biases should dominate. 



Our results are more similar to those of iNissen fc Schuster] (l2010l . l201ll ): ISchuster et al. 



(I2OI2I ). who find roughly equal-size populations distinct in age, abundances and orbits when 
they divide their (metal-rich, [Fe/H]> —1.6) stellar halo sample between low-a and high-a 
stars by drawing a line in the [Mg/Fe] vs. [Fe/H] plane. They interpret the properties of 
the low-a population (in which the stars are younger and on more eccentric orbits) as being 
consistent with an accreted population — and hence have an accreted fraction for the stellar 
halo at these metallicities similar to our own estimates. Note that we also find stars w' i th sig- 
nificantly lower [a/Fe] (less than zero) than any of the stars in the lNissen fc Schusterl (120101 ) 
sample, possibly because our large survey volume (out to 9 kpc from the Sun) encompasses 
debris from peculiar accretion events not represented in their local sample (within 335 pc of 
the Sun). 
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5.3. Conclusion 



We conclude that, while there are strong indications of distinct populations with distinct 
formation mechanisms in the stellar halo, definitiv e evidence for a significan t contribution 
from in situ/ kicked- out- disk stars remains elusive. iNissen &: Schusterl (120111 ) have demon- 
strated that using abundance patterns along with orbital distributions can play a key role in 
making this identification. Our own work points to the importance of mapping a significant 
volume of the halo in order to confirm that such local studies are representative of global 
properties. In particular, the volume of our survey has allowed us for the first time to find 
candidates for a metal-rich kicked- out- disk population. Quantifying the size of the kicked- 
out- disk population more generally could provide valuable constraints on the hydrodynamical 
models of galaxy formation. 

In addition, it is apparent that our selection of M giants with unusual, halo-like RVs picks 
out a stellar halo population dominated by accreted stars. Hence, searching for accretion 
events in RV-outlying samples of M giants should be particularly fruitful and motivates 
further interest in the putative groupings of stars found in our RV survey. In a companion 
paper we go on to explore wh at more the kinemat ic properties of these groupings could be 
telling us about their origins (I Johnston et al.ll2012[ ). 



Many thanks to Eric Bell, Paul Harding and Tim Beers for elucidating discussions at 
the January 2012 meeting of the American astronomical Society, which added a healthy dose 
of skepticism to the discussion of our results. The work of KVJ and AAS on this project 
was funded in part by NSF grants AST-0806558 and AST-1107373. AAS acknowledges 
support from The Vassar College Committee on Research and the Columbia Science Fellows 
Program. 
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Fig. 1. — Schematic showing the location of the three populations - in-situ, kicked- out- disk, 
and accreted - in the [a/Fe]-[Fe/H] plane. 
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Fig. 2. — Aitoff projection of the Galactic coordinates of the 1799 M giants observed at 
FMO and CTIO. The program began in February, 2005 and is ongoing. Stars were selected 
to primarily sample a wide sweep of the Milky Way's thick disk and nearby halo. 
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Fig. 3. — A typical medium-resolution spectrum for the program stars, with prominent lines 
labeled. The dominant features in M giants in the spectral band we are studying are the Mg 
b triplet, at 5167, 5173, and 5184 A, and the Na D doublet at 5889 and 5896 A. A number 
of strong Fe, Cr, and Ti lines/blends are also present. 
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Fig. 4. — Medium-resolution spectra of four stars observed at Fan Mountain Observatory, 
spanning a range of (J — ^^5)0 colors for typical program stars. The star IDs and colors are 
printed at the left. The topmost spectrum is of the type M4 giant HD101153. The increasing 
prominence of molecular bands is apparent for the redder (cooler) stars. 
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Fig. 5. — Radial velocities as a function of Galactic longitude. Panel (a) shows Uhei vs. / 
for all 1799 program M giants. Panel (b) shows Vqsr vs. / for the program stars. In panel 
(c), fGSR has been divided by cos(6) to accentuate differences between stars with disk-like 
and non-disk-like motions. Panel (d) repeats panel (c), where stars with high- resolution 
spectra obtained are overplotted in color. The majority of stars observed fall along the 
main disk trend, but several groupings of stars that deviate from this trend are conspicuous. 
The different colors used to mark the stars with high-resolution spectra correspond to these 
apparent groupings. 
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Fig. 6. — Teff-log g isochrones for populations of 3 Gyr (dotted line), 5 Gyr (solid line), and 
10 Gyr (dashed) populations. The range of log(Te//) derived for the APO/KPNO program 
stars is 3.556-3.602. At higher Tefr, the log g values begin to diverge significanty. 
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Fig. 7. — [Ti/Fe] vs. [Fe/H] for program stars with high-resolution spectra and S[Ti/Fe] <0.20 
dex, including the four red giant standards Arcturus, a Tau, u Vir, and 6 Vir (standard 
stars shown as yellow circles). Filled stars are color-coded according to possible groupings 
identified via radial velocities, consistent with the scheme used in Figure [51 The high- 
resolution M giants have been outlined as follows: bold triangles correspond to accreted stars; 
bold squares correspond to kicked- out- disk stars; bold diamonds could be accreted, kicked- 
out-disk, or in-situ-halo stars. Milky Way disk and field stars taken from the literature are 
shown in black and are coded as follows: triangles - Fulbright (2002); plus signs - Reddy et 
al. (2003); squares - Brewer & Carney (2006). For comparison with the red giant standard 
stars, five randomly selected thick disk stars - based on RVs from Fig. [5] - are shown as 
yellow stars. Typical error bars are shown in the lower left in black. 
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Fig. 8. — The heliocentric radial velocity as a function of height above the Galactic plane 
for the 34 program stars with high-resolution spectra. Same color-coding as used as in Fig. 
[5l The dashed lines show 2.5adisk and the dotted line indicates the scale height of the thick 
disk. 
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Fig. 9. — The heliocentric radial velocity as a function of Galactic longitude for: (a) program 
stars with b < —30° and (b) all program stars. The box from 30° < / < 7 5° indicates the 



location and radial velocity range of the Aquarius stream detected in the IWilliams et al. 
(I2OIII ) survey. 



Table 1. Observed Properties of Program Stars. 



ID 


I 


b 


{Ks)o 


{J - Ks)o 


Date 


Inst 


S/N Uhel,h/l'hel,m Grp 


o 




deg 


deg 






UT 




km 





1521021+082320 


12.1085 


49.9949 


7.458 


0.969 


1546044+061554 


14.3431 


43.6135 


7.745 


0.961 


1535178+135331 


22.2041 


49.6287 


6.569 


1.028 


1640358+060251 


22.6837 


31.7562 


5.655 


1.013 


1509307+252912 


37.8279 


59.1097 


8.680 


0.991 


1546468+270052 


43.1496 


51.1734 


6.965 


0.972 


1530257+323409 


51.7506 


55.2995 


5.997 


1.053 


1719596+301146 


53.2262 


31.8831 


6.863 


0.919 


1731343+401543 


65.3999 


31.7765 


5.335 


1.019 


0913092+450916 


175.2015 


43.4028 


8.045 


0.937 


0807206+435418 


176.0975 


31.6497 


6.552 


0.937 


0935556+414046 


179.7566 


47.7631 


6.179 


0.920 


0903471+414944 


179.7632 


41.7691 


6.322 


0.955 


0830244+283812 


194.9298 


33.0562 


6.585 


0.973 


1011498+230504 


210.2972 


53.7906 


7.908 


0.934 


0918573+145749 


215.4937 


39.3263 


8.780 


0.908 


0938204+112950 


222.2501 


42.1710 


8.349 


0.980 


1001593+114507 


225.6258 


47.4509 


8.399 


0.932 


0932038+055521 


227.7177 


38.1643 


7.397 


1.073 


0950147+082356 


227.8484 


43.3082 


4.981 


1.154 


1101590+084043 


243.2806 


58.2289 


8.748 


0.931 


1024571+004900 


243.5677 


46.1144 


7.173 


1.054 


1101024^003004 


254.5704 


51.6925 


7.208 


0.923 


1115376+000800 


258.5365 


54.5283 


8.248 


1.037 


1054035-065124 


258.7632 


45.7090 


7.355 


1.024 


1037414-121048 


259.0329 


39.0328 


5.977 


1.142 


1136527+025949 


263.2681 


59.9961 


7.387 


0.988 


1145072+013727 


268.2122 


59.9476 


7.274 


1.029 


1105038-164000 


269.3479 


39.1698 


7.386 


0.980 


1131243-055825 


269.6677 


51.6539 


6.526 


1.158 


1206183-035045 


281.8293 


57.1651 


6.521 


1.095 


1313176-164220 


310.4547 


45.8468 


8.848 


0.942 


1314125-132352 


311.4065 


49.0990 


7.125 


0.943 



Ol-Mar-2010 


ECHLR 


62 


-192.9/-183.3 


A 


a 


26-Feb-2010 


ECHLR 


75 


-268.6/-270.6 


A 


a 


30-Mar-2009 


ARCES 


114 


53.9/61.1 


d 


is/ko/a 


30-Mar-2009 


ARCES 


80 


35.4/49.4 


d 


is/ko/a 


25-Feb-2010 


ECHLR 


84 


-231.6/-231.6 


A 


a 


26-Fcb-2010 


ECHLR 


98 


-217.3/-219.6 


A 


a 


30-Mar-2009 


ARCES 


138 


-8.6/-4.6 


d 


is/ko/a 


02-Mar-2010 


ECHLR 


52 


-30.1/-27.0 


d 


is/ko/a 


26-Feb-2010 


ECHLR 


111 


-28.2/-27.8 


d 


is/ko/a 


26-Feb-2010 


ECHLR 


105 


-107.0/-108.2 


B 


a 


26-Feb-2010 


ECHLR 


123 


-122.5/-131.7 


B 


ko 


26-Feb-2010 


ECHLR 


174 


-252.7/-248.3 


B 


a 


26-Fcb-2010 


ECHLR 


195 


-204.8/-210.6 


B 


a 


25-Feb-2010 


ECHLR 


174 


264.5/ 269.9 


C 


a 


02-Mar-2010 


ECHLR 


89 


238.4/242.4 


C 


a 


02-Mar-2010 


ECHLR 


79 


229.9/235.2 


C 


is/ko/a 


25-Feb-2010 


ECHLR 


97 


264.2/265.6 


C 


a 


02-Mar-2010 


ECHLR 


77 


179.2/183.2 


C 


is/ko/a 


02-Mar-2010 


ECHLR 


82 


190.3/171.7 


C 


ko 


02-Mar-2010 


ECHLR 


107 


143.8/167.4 


c 


ko 


02-Mar-2010 


ECHLR 


73 


263.8/263.9 


D 


a 


25-Feb-2010 


ECHLR 


109 


304.5/314.1 


D 


a 


25-Feb-2010 


ECHLR 


172 


264.9/267.0 


D 


is/ko/a 


02-Mar-2010 


ECHLR 


72 


196.4/194.6 


D 


ko 


25-Fcb-2010 


ECHLR 


154 


267.1/269.9 


D 


is/ko/a 


25-Fcb-2010 


ECHLR 


185 


248.2/257.1 


D 


ko 


25-Fcb-2010 


ECHLR 


117 


162.1/159.3 


D 


a 


26-Feb-2010 


ECHLR 


105 


146.6/145.9 


D 


a 


02-Mar-2010 


ECHLR 


89 


349.7/354.6 


D 


is/ko/a 


26-Feb-2010 


ECHLR 


129 


123.2/137.9 


D 


is/ko/a 


25-Feb-2010 


ECHLR 


151 


220.6/238.9 


D 


a 


Ol-Apr-2009 


ARCES 


82 


113.7/96.9 


E 


ko 


Ol-Apr-2009 


ARCES 


132 


99.3/98.9 


E 


ko 



Table 1 — Continued 



ID 


I 


b 


{Ks)o 


(J - Ks)o 


Date 


Inst 


S/N 


«hel,h/fhel,m 


Grp 


o 




deg 


deg 






UT 






km s~^ 






1334567-055158 


322.2190 


55.3680 


6.968 


1.035 


Ol-Apr-2009 


ARCES 


100 


71.5/77.6 


E 


ko 


1355160-053350 


330.5657 


53.8480 


7.352 


1.057 


02-Mar-2010 


ECHLR 


86 


91.9/96.2 


E 


is/ko/a 


1442035-162350 


337.8142 


38.8681 


7.298 


1.050 


02-Mar-2010 


ECHLR 


73 


152.2/145.9 


E 


is/ko/a 


1404515-004157 


338.2993 


57.0434 


6.590 


1.109 


02-Mar-2010 


ECHLR 


73 


114.9/121.4 


E 


is/ko/a 


1431302-051730 


343.4025 


49.5528 


8.156 


0.909 


Ol-Apr-2009 


ARCES 


62 


162.0/156.5 


E 


a 


1518065-115536 


350.0932 


37.1732 


7.588 


1.045 


30-Mar-2009 


ARCES 


91 


185.2/183.0 


E 


is/ko/a 



Table 2. Equivalent Width Measurements. 



ID 7443.02 7447.38 7461.52 7498.53 7507.26 7511.02 7531.14 7540.43 7547.91 7568.89 7583.79 7489.57 7496.12 7474.94 

AAAAAAAAAAAAAA 

Fel Fel Fel Pel Fel Fel Pel Fel Fel Pel Fel Til Til Til 



1521021+082320 


53.5 


34.5 


83.1 






162.4 


100.1 


56.9 


20.0 


94.2 


147.0 


56.5 


42.1 


26.7 


1546044+061554 


50.0 


34.2 


86.1 


33.8 


66.2 


171.6 


98.2 


51.0 




92.0 


144.4 


48.6 


33.5 


27.9 


1535178+135331 


49.1 


42.8 


104.2 


42.8 




165.8 


109.3 


46.2 


22.7 


99.0 


156.6 


104.7 




30.0 


1640358+060251 


77.5 


51.9 


108.5 


52.5 


95.6 


186.9 


119.4 


67.5 


31.0 


108.4 


159.8 


113.8 


87.5 


33.0 


1509307+252912 


40.1 


31.2 


80.4 


28.1 






94.2 


43.4 


17.8 




130.8 


35.5 




30.0 


1546468+270052 


64.0 


36.8 


98.5 


45.4 


87.2 


180.6 


115.3 


58.1 


31.5 


108.9 


157.5 


72.4 


49.3 


21.1 


1530257+323409 


63.0 


47.1 


103.8 


52.6 


70.6 


171.1 


101.3 


60.0 


24.2 


94.6 


147.0 


115.0 


102.5 


32.3 


1719596+301146 


71.3 


57.0 


85.3 




89.8 


169.5 


114.9 


65.0 


31.8 


109.5 


155.2 


108.7 


78.1 


41.7 


1731343+401543 


81.3 


60.4 


120.6 


64.1 


108.5 


201.0 


127.3 


76.0 


35.9 


123.2 


175.6 


138.2 


107.8 


50.6 


0913092+450916 


62.9 


41.6 


90.3 


40.9 


77.0 


166.5 


108.0 


47.9 




99.8 


143.9 


60.8 


37.9 


24.4 


0807206+435418 




68.6 


124.4 


47.7 


118.6 


205.3 


134.9 


89.1 


41.4 


134.2 


181.6 


139.0 


100.4 


56.9 


0935556+414046 


46.2 


26.3 


67.7 


27.3 


59.4 


154.2 


88.3 


45.0 


19.9 


82.2 


130.8 


35.9 


22.1 


13.0 


0903471+414944 


60.0 


43.5 


91.8 


54.3 


77.0 


161.1 


102.0 


47.3 


20.0 


96.4 


140.8 


88.2 


64.8 


Sij30.3 


0830244+283812 


49.0 


34.5 


85.0 


32.9 


66.3 


159.3 


100.1 


46.8 


14.2 


91.8 




61.0 


47.2 


1 16.9 


0918573+145749 


51.9 


21.1 


76.1 




64.0 


157.3 


94.7 


36.7 


12.6 


88.3 




39.3 


24.7 


19.0 


0938204+112950 


60.5 


37.4 


94.6 


39.1 


78.5 


165.4 


107.3 


54.4 


21.6 


105.4 




72.8 


61.9 


20.2 


1001593+114507 


66.8 


46.4 


94.3 




87.0 


182.6 


109.5 


55.0 


17.4 


97.3 


150.9 


87.8 


69.4 


43.1 


0932038+055521 


67.5 


45.0 


94.1 


43.8 


74.7 


162.9 


102.3 


59.7 


33.8 


98.7 


149.0 


105.7 


97.8 


38.9 


0950147+082356 


55.6 


45.1 


101.7 


47.3 




156.7 




55.9 


22.6 


90.2 


140.3 


108.5 


96.1 


26.1 


1011498+230504 


63.3 


44.7 


95.7 


52.4 


77.6 


171.6 


105.5 




27.6 


95.0 


152.1 


67.0 


45.8 


30.7 


1101590+084043 


41.9 




71.8 




60.0 


143.0 


100.7 


44.2 


18.7 


80.6 


135.7 


43.8 


33.0 


16.0 


1024571+004900 


58.0 


37.3 


105.8 


41.6 


76.9 


169.5 


114.8 


60.0 


21.7 


97.6 


158.1 


82.8 


65.9 


26.9 


1101024-003004 


42.2 


25.1 


81.8 


24.8 


61.3 


155.4 


93.8 






82.4 


136.3 


43.3 


30.1 


12.1 


1115376+000800 


63.6 


45.5 


95.3 




74.0 


162.8 




60.2 


26.5 


96.4 


145.5 


104.7 


89.0 


27.8 


1054035-065124 


59.1 


35.4 


97.0 


35.3 


74.0 


158.0 


100.6 




17.5 


89.1 


147.6 


100.8 


79.1 


22.8 


1037414-121048 


65.3 




102.7 


39.6 


70.3 


151.8 


92.5 


54.4 


21.2 


94.2 




103.1 


87.1 


29.2 


1136527+025949 


44.5 


26.3 


81.3 


27.4 


66.2 


162.8 


97.3 


39.9 


15.6 


85.6 


146.6 


47.0 


35.0 


17.0 


1145072+013727 


55.6 




99.9 




85.2 


171.2 


110.0 


60.9 


25.5 


93.4 


145.5 


90.7 


74.3 




1105038-164000 


45.8 


20.9 


81.1 


21.7 


76.2 


170.9 


106.6 


37.9 




87.3 


150.1 


65.8 


43.0 


17.7 


1131243-055825 


45.2 


28.2 




29.1 


67.3 


134.5 




50.2 


21.9 


59.8 


132.1 


117.9 


87.4 


36.9 


1206183-035045 


49.6 


36.9 


94.3 


42.5 


72.4 


158.5 


101.6 


56.5 


19.6 


99.7 


152.4 


91.7 


72.3 


22.9 


1313176-164220 


66.9 


50.1 


103.6 


53.5 


87.3 


180.6 


115.5 


66.3 


43.8 


98.0 


156.8 


110.0 


96.8 


24.8 



Table 2 — Continued 



ID 


7443.02 


7447.38 


7461.52 


74QS 




7511.02 


7531.14 


7540.43 


7547.91 


1 ouo.oy 


7^8^ 7Q 


7489.57 


7496.12 


7474.94 




X 


A 


A 


X 


X 

l\ 


A 


A 


X 


X 


A 


A 


A 


X 


A 






X'tJi 




IT til 




X'tJi 


X'tJi 




I? til 






Til 


Til 


Til 


1314125-132352 


59.6 


43.0 


93.0 


52.1 


84.1 


173.5 


106.4 


56.5 


26.7 


97.2 


144.0 


89.5 


70.4 


22.6 


1334567-055158 


65.8 


50.7 


95.7 


48.5 


79.9 


170.6 


110.7 


67.2 


25.2 


97.1 


149.7 


107.8 


82.7 


22.0 


1355160-053350 


50.0 


32.0 


73.1 


25.1 


57.0 


146.5 


89.3 






80.6 


149.4 


69.0 


64.4 


53.8 


1442035-162350 


42.5 


33.4 


87.6 


32.4 


64.9 


154.1 


91.9 


46.2 


13.3 


87.3 


133.6 


74.6 


60.9 


40.5 


1404515-004157 


56.2 


33.2 


112.4 


45.1 


70.6 


171.2 


113.5 


63.2 


20.6 


104.3 


167.4 


107.7 


87.6 


26.2 


1431302-051730 


45.9 


28.5 


85.4 


31.3 


70.3 


164.4 


104.6 


39.5 




86.2 


140.3 


39.2 


28.0 




1518065-115536 


54.1 


43.9 


97.6 


43.4 


77.4 


162.0 


97.3 


56.2 


17.1 


91.5 


141.9 


84.3 


68.4 


20.2 


Arcturus 


65.5 


45.5 


94.6 


31.0 


89.3 


185.0 


115.8 


48.0 


31.6 


106.6 


148.4 


54.1 


72.9 


35.9 


a Tau (KPNO nl) 


93.0 


63.4 


120.8 


78.4 


115.4 


212.7 


132.4 


82.5 


45.3 


129.0 


179.0 


127.9 


109.1 


99.6 


a Tau (KPNO n4) 


95.7 


65.8 


120.0 


73.5 


113.9 


211.0 


132.2 


78.5 


45.8 


121.6 


181.0 


124.4 


104.0 


90.4 


5 Vir (APO) 


79.1 


57.4 


115.9 


65.7 


89.4 


193.3 


116.0 


73.9 


35.6 


112.2 


153.3 


126.8 


121.3 


99.3 


& Vir (KPNO) 


76.1 


54.3 


117.4 


68.0 


95.6 


184.0 




69.7 


37.9 


112.0 


164.2 


131.5 


118.8 


99.9 


1/ Vir (APO) 


73.4 


50.7 


101.7 


59.2 


85.7 


171.8 


110.7 


60.9 


29.3 


98.7 


155.4 


109.3 


93.4 


CO 76.3 
00 


V Vir (KPNO) 


66.9 


51.2 


97.0 


57.2 


82.8 


175.3 


107.4 


60.3 


25.4 


93.0 


149.9 


103.9 


92.5 


1 76.4 



Table 3. Derived Properties of Program Red Giants. 



ID 


T a 
-'eff 


log g 


c 






pfFp/Hl 

CI J. C/ iJ-J 




[Ti /FpI 


pfTi /FpI 

CI -L 1/ J. CI 




1521021+082320 


3900 


0.4 


1.48 


6.50 


-0.95 


0.08 


3.67 


-0.28 


0.07 


4.7 


1546044+061554 


3900 


0.5 


1.58 


6A4 


-1.01 


0.11 


3.76 


-0.13 


0.17 


5.6 


1535178+135331 


3850 


0.3 


1.62 


7.03 


-0.42 


0.17 








2.7 


1640358+060251 


3850 


0.9 


1.57 


7.24 


-0.21 


0.16 


4.51 


-0.18 


0.05 


1.6 


1509307+252912 


3900 


0.4 


1.28 


6.35 


-1.10 


0.13 








8.6 


1546468+270052 


3900 


0.7 


1.64 


6.83 


-0.62 


0.09 


4.09 


-0.19 


0.15 


3.3 


1530257+323409 


3800 


0.5 


1.39 


6.84 


-0.61 


0.16 


4.68 


0.39 


0.12 


2.5 


1719596+301146 


3900 


1.0 


1.35 


7.16 


-0.29 


0.13 


4.71 


0.10 


0.12 


2.3 


1731343+401543 


3800 


1.1 


1.71 


7.45 


0.00 


0.18 


4.89 


-0.01 


0.04 


1.0 


0913092+450916 


3950 


0.7 


1.48 


6.74 


-0.71 


0.10 


3.96 


-0.23 


0.08 


5.4 


0807206+435418 


3900 


1.3 


1.69 


7.52 


0.07 


0.12 


4.91 


-0.06 


0.07 


1.4 


0935556+414046 


4000 


0.5 


1.31 


6.37 


-1.08 


0.12 


3.73 


-0.09 


0.13 


2.5 


0903471+414944 


3900 


0.6 


1.42 


6.72 


-0.73 


0.16 


4.24 


0.07 


0.04 


2.7 


0830244+283812 


3900 


0.4 


1.54 


6.41 


-1.04 


0.10 


3.88 


0.02 


0.15 


3.2 


1011498+230504 


4000 


0.7 


1.63 


6.57 


-0.88 


0.08 


4.10 


0.08 


0.12 


5.0 


0918573+145749 


4000 


0.4 


1.75 


6.19 


-1.26 


0.12 


3.87 


0.23 


0.13 


9.0 


0938204+112950 


3900 


0.5 


1.57 


6.63 


-0.82 


0.08 


4.11 


0.03 


0.13 


7.1 


1001593+114507 


4000 


0.5 


1.74 


6.51 


-0.94 


0.13 


4.45 


0.49 


0.18 


6.6 


0932038+055521 


3700 


0.9 


1.20 


7.33 


-0.12 


0.11 


4.73 


-0.05 


0.13 


3.2 


0950147+082356 


3600 


0.5 


1.37 


7.12 


-0.33 


0.17 


4.48 


-0.09 


0.16 


1.7 


1101590+084043 


4000 


0.4 


1.47 


6.31 


-1.14 


0.13 


3.82 


0.06 


0.16 


8.1 


1024571+004900 


3800 


0.3 


1.57 


6.59 


-0.86 


0.08 


4.07 


0.03 


0.12 


4.3 


1101024-003004 


4000 


0.5 


1.58 


6.22 


-1.23 


0.10 


3.82 


0.15 


0.12 


4.2 


1115376+000800 


3800 


0.7 


1.33 


6.95 


-0.50 


0.08 


4.57 


0.17 


0.11 


6.5 


1054035-065124 


3800 


0.0 


1.42 


6.46 


-0.99 


0.13 


4.25 


0.34 


0.10 


4.9 


1037414-121048 


3650 


0.1 


1.29 


6.83 


-0.62 


0.14 


4.41 


0.13 


0.13 


3.3 


1136527+025949 


3950 


0.4 


1.69 


6.27 


-1.18 


0.10 


3.79 


0.07 


0.14 


4.8 


1145072+013727 


3800 


0.4 


1.51 


6.72 


-0.73 


0.08 


4.26 


0.09 


0.17 


3.4 


1105038-164000 


3900 


0.4 


2.04 


6.13 


-1.32 


0.08 


3.87 


0.29 


0.15 


4.8 


1131243-055825 


3650 


0.0 


1.23 


6.52 


-0.93 


0.19 


4.40 


0.43 


0.11 


3.7 


1206183-035045 


3700 


0.3 


1.47 


6.61 


-0.84 


0.10 


4.07 


0.01 


0.12 


3.7 


1313176-164220 


3900 


0.9 


1.55 


7.07 


-0.38 


0.15 


4.64 


0.12 


0.10 


6.2 


1314125-132352 


3950 


0.8 


1.41 


6.89 


-0.56 


0.11 


4.49 


0.15 


0.07 


3.1 


1334567-055158 


3750 


0.8 


1.36 


7.16 


-0.29 


0.09 


4.69 


0.08 


0.08 


3.1 


1355160-053350 


3800 


0.0 


1.47 


6.26 


-1.19 


0.16 


3.88 


0.17 


0.19 


5.1 



Table 3 — Continued 



ID 


T'efl 


log 9 




A(Fe) 


[Fe/H] 


e[Fe/Hl 


A(Ti) 


[Ti/Fe] 


e[Ti/Fe] 


dist (kpc) 


1442035-162350 


3800 


0.2 


1.50 


6.36 


-1.09 


0.12 


3.93 


0.12 


0.11 


4.8 


1404515-004157 


3750 


0.1 


1.69 


6.52 


-0.93 


0.11 


4.29 


0.32 


0.16 


3.6 


1431302-051730 


4000 


0.5 


1.71 


6.35 


-1.10 


0.10 


3.78 


-0.02 


0.25 


6.2 


1518065-115536 


3800 


0.3 


1.44 


6.57 


-0.88 


0.13 


4.17 


0.15 


0.11 


5.5 


Arcturus 


4250 


1.4 


1.76 


6.82 


-0.63 


0.09 


4.56 


0.29 


0.11 




a Tau(KPNO nl) 


3900 


1.3 


1.59 


7.50 


0.05 


0.11 


4.98 


0.03 


0.12 




a Tau(KPNO n4) 


3900 


1.3 


1.60 


7.62 


0.17 


0.10 


4.99 


-0.08 


0.06 




5 Vir{KPNO) 


3700 


0.9 


1.39 


7.63 


0.18 


0.13 


5.11 


0.07 


0.11 




5 Vir(APO) 


3700 


0.8 


1.35 


7.62 


0.17 


0.14 


5.05 


0.03 


0.15 




V Vir(KPNO) 


3750 


1.0 


1.40 


7.11 


-0.34 


0.15 


4.58 


0.02 


0.12 




V Vir(APO) 


3800 


0.8 


1.36 


7.19 


-0.26 


0.14 


4.64 


0.03 


0.03 
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Table 4. Standard Star Chemical Abundance Comparisons 



star [Fe/H] [Fe/H]s&L [Ti/Fe] [Ti/Fe]s&L 



Arcturus -0.63±0.09 -0.67±0.11 0.29±0.11 0.41±0.12 

a Tau O.llitO.lO 0.07 -0.04±0.08 -0.03 

u Vir -0.30±0.15 -0.02±0.20 0.01±0.20 0.06±0.17 

<5Vir 0.18±0.17 0.13±0.17 0.00±0.13 0.07±0.20 



